ON UNBIASED STOCHASTIC NAVIER-STOKES 

EQUATION 
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Abstract. A random perturbation of a deterministic Navier-Stokes 
equation is considered in the form of an SPDE with Wick type non- 
linearity. The nonlinear term of the perturbation can be characterized 
as the highest stochastic order approximation of the original nonlinear 
term itVu. This perturbation is unbiased in that the expectation of a so- 
lution of the perturbed equation solves the deterministic Navier-Stokes 
equation. The perturbed equation is solved in the space of generalized 
stochastic processes using the Cameron-Martin version of the Wiener 
chaos expansion. It is shown that the generalized solution is a Markov 
process and scales effectively by Catalan numbers. 



1. Introduction 
In this paper we will consider a deterministic Navier-Stokes equatiorQ 

(1.1) d t u (t,x) = di (a ij (t,x)djU {t,x)) 

- v% (i, x) d k u (t, x) + VP (t, x) + f (t, x) , 
uo (0, x) = w(x), div uo = 0, 

and its stochastic perturbations: 

(1.2) dtv (t, x) = di (a ij (t, x) djv (t, x)) - v k (t, x) d k v(t, x) 

+ VP (t,x) + i(t,x) 

+ [a% x)d,iy (t, x) + g (t, x) - VP (t, x)]W t , 
v (0, x) = w (x) , div v = 0, 
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and 

(1.3) d t u(t,x) = di (a ij (t,x)d j u(t,x)) - u k (t, x) <)d k u(t, x) 

+ VP (t, x) + f (t, x) 

[<7*(t,a;)^u(t,x) + g(t,x)- VP(t,i)]lf ( , 
u (0, x) = w (x) , div u = 0, 

where < t < T, x G R d , d > 2, and W* is a cylindrical Wiener process in a 
separable Hilbert space Y. The coefficients a lJ ,a l and the functions f, g are 
deterministic, a 1 and g are K-valued. Symbol stands for Wick product 
(see Section T2.2.1I and references [3J, [29j). Wick product is a stochastic 
convolution. It could be interpreted as a generalized Malliavin divergence 
operator with respect to Gaussian measure associated with white noise W 
(see m)- 

Stochastic PDEs involving Wick product type nonlinearity were originally 
discussed in the literature related to the Parisi-Wu program (see [26], [12] 
and also [20J (Section 6)). In these papers Wick product was defined by 
(Gaussian) invariant measures for the related PDEs. Other related papers 
include: |2J, [5], \T$, [E]-[I7], [H], etc.). 

Equations (jl.2p and (|1.3p are stochastic perturbations of the deterministic 
Navier-Stokes equation (jl.ip . It is shown in Section 3 that the generalized 
mean, i.e. the zero-order coefficient in the Wiener chaos expansion (|1 .8j) of 
the solution of equation (jl.3p . is a solution of equation (jl.ip . i.e. 

(1.4) Eu(t,x) =u (t,x), 

where uo (i, x) is a solution of equation (jl.ip . In other words, the solution 
of stochastic Navier-Stokes equation (jl.3p is a mean preserving (unbiased) 
random perturbation of deterministic Navier-Stokes equation (jl.ip . 

Obviously, this nice property does not hold for equation (jl.2p or other 
standard stochastic perturbations of Navier-Stokes equation (e.g. random 
initial conditions, random forcing, etc.) 

In fact, equation (jl.3p could be viewed as an approximation of stochastic 
Navier-Stokes equation (jl.2p . Indeed, under certain natural assumptions, 
the following equality holds: 



oo 



(1.5) vVv = E ""v0yW 

n=0 

where T> n is the n th power of Malliavin derivative T>. Taking into account 
expansion (jl.5p . v<}Vv could be viewed as an approximation of the product 
vVv. In fact, v<)Vv is the highest stochastic order approximation of vVv 
(see Appendix I, Proposition |4] and Remark II ip . 

Stochastic Navier-Stokes equation (jl.2p is reasonably well understood and 
there exists substantial literature on its analytical properties as well as its 
derivation from the first principles (see e.g. [23], [23] and the references 
therein). In this paper we will be focusing mostly on equation (jl.3p . 
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Burger's equation with Wick product was considered in [6], [8], and [9], 
see also the references therein. 

It was shown in [24] that under reasonable assumptions stochastic Navier- 
Stokes equation (jl.2j) has a square integrable solution. Moreover, this solu- 
tion can be formally written in the Wiener chaos expansion form: 

a 

where a G </} is the Cameron-Martin basis generated by Wt, v a (t, x) = 
E (v (t, x) £ a ), and J is the set of multiindices a = {ct k , k > 1} such that 
for every k, a k S No(No= {0, 1,2,.. .}) and \a\ = ^2 k a k < oo. 

It was shown in [23] that the Wiener chaos coefficients v a (t, x) satisfy 
the propagator equation: 

(1.6) d t v a (t, x) = di [a ij djv a (t, x)) - VP (t, x) + f (t, x) I {H=0} 

(t,x) 

V 0</3<a 

^2Vak[(<T\e k ) Y diV a ( k ) (t,x) + I { \ a \ =1} (g,e k ) Y ]; 

k 

v Q (0,x) = w Q (a;),div v a = 0, 
where a(k) = (ct%, a.2, ...a k _i,a k — l,a k +i, ...) and 



One advantage of the Wiener chaos representation is that it provides con- 
venient explicit formulae for computing statistical moments of the random 
field u(t,x) via Wiener chaos coefficients (see [23], [24]). For example, 

Eu i {t,x)=ui(t,x), E(^(t,x)n J (t,y)) = £ H<00 <(*, x)v? a (t, y), 

In this paper, we prove that the WCE coefficients of a solution of the 
unbiased stochastic Navier-Stokes equation (|1.3p are given by 

(1.7) d t u a (t, x) = di (a ij dju a (t, x)) - VP (t, x) + f (t, x) / {|q |=o} 

0</3<a V 

+ ^2 V^li^i 1 ' x )> e k{t)) Y d j u a(k) (t, X) 
k 

+ (g (t,x) ,e k {t)) Y 1 H=1 ] 
u Q (0,x) = w(x),div u a = 0. 
Clearly, this system of equations is much simpler than equation (jl.6p . 
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If a = 0, then u a (t,x) is a solution of deterministic Navier-Stokes equa- 
tion (jl.ip . The remaining components are governed by Stokes equations and 
could be solved sequentially. From the computational point of view this is 
a substantial advantage. Indeed, the propagator for equation (|1.2j) is a full 
nonlinear system while equation (|1.3p is a lower triangular system and only 
the first equation of this system is nonlinear. 

An important feature of equation fjl .3f) is that 

(u k (t)0d k u(t),u(t)) 

Therefore, one could not expect a solution of (11.30 to be square integrable. 
This effect is not specific to stochastic Navier-Stokes equation. In fact, it is 
common for a large class of stochastic bilinear PDEs (see e.g. [15] . [16]). 

In this paper we consider equation (jl.3j) in the class of formal Wiener 
chaos expansions and show that a formal series 

(1.8) u(t,x) = y^Ug (t,x)£ a 

solves (II. 3p if and only if u a (t, x) are given by equation ([1.7]) . To make this 
solution square integrable, we rescale it using second quantization operators 
(see Appendix I, 5.1). It is shown that u(t, x) is the limit of square integrable 
solutions of the rescaled equations (see Proposition [2}. 

Convergence of this solution is determined by a system of positive weights 

{ r «}|a|<oo such that 

(1-9) \\An-= Yl r «ll Uc *(*)llL((0,T);R<*) < 00 ■ 

\a\ <oo 

It turned out, that Catalan numbers (see [30] . [9]) are critical for an appro- 
priate choice of the weights r a in fjl .9]) (see Proposition [I]). 

In addition, it was shown that a solution of equation (jl.3p belongs to the 
intersection of Sobolev spaces H2 (R d ) n Ml (R d ) for p > d. We have also 
demonstrated that uniqueness of a solution of equation (j 1 . 3 1) holds under 
the same assumptions that guarantee uniqueness for the related determin- 
istic Navier-Stokes equation. Although £ a in (|1.8p are not (J 7 t H/ )-adapted, 
we prove that the generalized solution is (J 7 i w/ )-adapted and Markov (see 
Theorem [2] and Corollary [5]) . 

It is not clear how, if at all, the unbiased Navier-Stokes equation fits 
into classical fluid mechanics. Nevertheless, equation (II. 3p is "physical" 
in that it could be derived from the second Newton law (under appropri- 
ate assumptions on the velocity field), much the same way as the classical 
Navier-Stokes equation (see Appendix I, 5.2). Also, it was shown recently 
(see [11] ) that, after Catalan type rescaling, finite dimensional projections 
of unbiased Navier-Stokes equation present an accurate and numerically in- 
expensive approximation of stochastic Navier-Stokes equation (II. 2p 
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We conclude this section with an outline of some notations that will be 
used in the paper. 

1.1. Notation. Let us fix a separable Hilbert space Y. The scalar product 
of x, y G Y will be denoted by (x, y)y- 

If u is a function on Ti d , the following notational conventions will be used 
for its partial derivatives: diU = du/dxi,dfj = d 2 u/dxidxj, dfU = du/dt, 
and Vu = du = {d\u, . . . , dga), and d 2 u = (dfjic) denotes the Hessian 
matrix of second derivatives. Let a = (ai, ay) be a multi-index, a% G 
N = {0, 1, 2, . . .} , i = 1, . . . , d, then d« = U d =1 d%. 

Vector fields on R d are denoted by boldface letters. This convention also 
applies if the entries of the vector field are taking values in a Hilbert space. 

We denote N = {1,2,...}. 

For a Banach space E, we denote C (\0,T],E) the space of continuous 
E- valued functions. 

Cg° = C^°(R d ) denotes the set of all infinitely differentiable functions on 
H d with compact support. 

For s G (-00,00), write A s = A s x = ( 1 - Ya=i d /dx 2 J . For p G 

[1, 00) and s G (—00, 00), we define the space = Hp(R d ) as the space of 
generalized real valued functions u with the finite norm 

|^|s,p — I A 

where | • \ p is the L p norm. Obviously, = L p . Note that if s > is an 
integer, the space Hp coincides with the Sobolev space W p = W p (R d ). 

The spaces Co°(R c! ), (R d ) can be extended to vector functions (de- 
noted by bold-faced letters). For example, the space of all vector functions 
u = (u 1 , . . . , u d ) such that A s u l G L p , I = 1, . . . , d, with the finite norm 

|uk P = (Ertp) 1/p ' 

is denoted by U s p = U s p (R d ). Similarly, we denote by M s p (Y) = M s p (H d ,Y) 
the space of all vector functions g = (g l )i<i<d, with l"-valued components 
g l , 1 < I < d, so that ||g|| s ,p = (Y^iWlstf) 1 ^ < 00 • Also, for brevity, the 
norm ||g||o,p is denoted by ||g|| p . 

When s = 0, W p (Y) = L P (Y) = 7L p (R d , Y). To forcefully distinguish 
Lp-norms in spaces of Y-valued functions, we write || • while in all other 
cases a norm is denoted by |-| . The duality (-, -) s between W s q (R d ) , and 
M.~ s (R d ) where p > 2 and q = pj (p — 1) is defined by 

d „ 

(<t>,tl>). = (<t>,il>)s P = Y l (A s ^)(x)(A-^)(x)dx,0G M^,VGIH- S . 
i=i J * d 
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2. Generalized Random Variables and Processes 

2.1. Wiener Chaos. To begin with, we shall introduce some basic notation 
and recall a few fundamental facts of infinite-dimensional stochastic calculus. 
Let us fix a separable Hilbert spaces Y and H = L2QO, T], y). Let {£i,i > 
1} be a complete orthonormal basis (CONS) in Y and{mj,i > 1} be a 
CONS in L 2 ([0,T]) . Denote by B the class of all CONS in H of the form 
{&k = efe(s) = mk 1 (s)£k 2 } and such that for each k, sup 0<s<T \mk(s)\ < 00. 
Obviously, for each k, sup < s <r |e^(s)|y < 00. Let us fix a CONS b = 
{e k ,k>l}£B. 

Let (CI, J- w , P) be a probability space with a cylindrical Brownian motion 
Wt in Y and J- w be the cr-algebra generated by W . Let ¥ w be the right 
continuous filtration of cr-algebras (F^)t>o generated by Wt. All the a- 
algebras are assumed to be P-completed. Hence 

00 

fc=i 

where \w\, k > l} is a sequence of independent standard one-dimensional 
Brownian motions in (Cl,^, P). We write W(ek) = L ek(t)dWt- For = 
e k (s) = m kl (s)£ k2 , 

W(e k )= [ e k (t)dW t = [ m kl (t)dw k2 (t), 
Jo Jo 

and 

00 „ t 

(2.1) W t = Y j ( e k (s)ds)W(e k ),0<t<T. 

k=i Jo 

Let a = {a k , k > 1} be a multiindex, i.e. for every k, a k £ No = 
{0, 1,2,.. .}. We shall consider only such a that \a\ = Yl k a k < 00, i.e., only 
a finite number of a k is non-zero, and we denote by J the set of all such 
multiindices. For a, (3 E J, we define 

a + P = (at +Pi,a 2 +P2,...), a\ = JJ a k \. 

k>l 

By e k we denote the multi- index a with a k = 1 and aj = for j ^ k. Write 

(2.2) a(k)=a-e k 

For a G J, write ff a := n^Li -^^(^ (ejfc)), where i/ n is the n* ft Hermite 
polynomial defined by H n (x) = (-l) N [d n e' x2 / 2 / dx n ^j e* 2 ' 2 . 
Let £ a = H a / \fa\. 

Theorem 1 (Cameron and Martin [lj). The set 5 = {£ Q = ^ Q (b), a G J} is 
an orthonormal basis in L 2 (CI, F w , P) , where T w is the a -algebra generated 
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by W. If E is a Hilbert space, 77 G L2 J 1 " 1 ^, P;-E) and r\ a = E(r/£ Q ), i/ien 

The expansion 77 = X^og j n at,a is often referred to as Wiener chaos ex- 
pansion. 

Remark 1. The basis £ a , a G J, can be obtained by differentiating stochastic 
exponent. Let Z be the set of all real-valued sequences z = (zk) such that 
only finite number of Zk is not zero. For a G J, denote d z = Hkd ak / {dzk)°' k 
and let 

e z = e z {t) = ^z fc e fc (t),0 <t<T, 

k 

(2.3) pt(z) = p t (e z ) = pt(z,h) = exp I J e z (s)dW s - ^ J \e z {s)\yds 
p( z ) = PT ( z ),z G Z,0 < t < T. 

It is a standard fact (see, for example, that H a = d^p(z)\ z= o,^ 

H a /y/a\. Since p(z) is analytic, it follows by H^W), 



(2.4) p{z) = p{z, b) = £ = £ -7=7^ 

a a v 

2.2. Generalized random variables and processes. Let b EB, and £ a = 

£ a (b), a G J. Let 

V = V(b) 

= fv = ^2 v aia ■ v a G R and only finite number of v a are not zero | . 

Definition 1. A generalized D-random variable with values in a convex 
topological vector (linear) space E with Borel a-algebra is a formal series 
u = ^2 a u a ^ a , where u a G E, £ Q = £ a (b), and h = {ek,k > 1} G B is a 
CONS in H = L 2 ([0, T],Y). 

Denote the vector space of all generalized P-random variables by V = 
T>'(h) = D'(h;E). The elements of T> are the test random variables for T> ! . 
We define the action of a generalized random variable u on the test random 
variable v by (u,v) = ^2 a v a u a . 

For a sequence u n G T>' and u G T>', we say that u n — > u, if for every 
v £ V, (u, v n ) — > (u, v) . This implies that u n = ^2 a u^ a — > u = Y^ a u a £, a if 
and only if u™ — > u a as n — > 00 for all a. 

Remark 2. Obviously, if u = ^2 a u a ^ a G V'(h;E), F is a vector space and 
f : E — )• F is a linear map, then 

f{u) = Y,f{u a )i a £V'{h-F). 
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Definition 2. An E -valued generalized D- process u(t) in[0,T] is aT>'(h; E)- 
valued function on [0, T] such that for each t G [0, T] 

u(t) = Y,Mt)ta£V'{b;E); 

a 

and u a (t) are deterministic measurable E -valued functions on [0, T]. We 
denote the linear space of all such processes by V(b;[Q, T], E). If E is a 
topological vector space and a generalized T>-process u(t) is continuous we 
write u G CT>'([0, T], b, E) (note that u(t) is continuous if and only if all 
coefficient functions u a are continuous in E. 

If there is no room for confusion, we will often say D-process (D-random 
variable) instead of generalized V -process (generalized V -random variable). 

If E is a normed vector space, we denote 
L^V'ih-^TlE)) 

= {«(*) = Y] «a(t)&, G £>'(*>; [0, T],E) : f \u a (t)\ E dt < oo, a G J}. 

Jo 

For u(t) = Y. a u a (t)£ a G Li(V(b; [0,T],E)) we define f*u(s)ds,0< t < T, 
■mV'(b;[0,T],E) by 



y u(s)ds = ^2 (^J u a ( 



s)ds ) f a ,0 < t < T. 



If = Ea«a(*)^a G ^'(b; [0,T],.E), then is differentiable in t if 
and only if u a {t) are differentiable in t . In that case, 

d_ 

It' 



: u(t) = u(i) = J^UaWta G V'([0,T],h,E). 



Example 1. A cylindrical Wiener process Wt,0 < t < T, in a Hilbert 
space Y, and its derivative dWt/dt = Wt are generalized Y -valued stochastic 
processes. Indeed, by \2. 

rt 

e k (s)ds£ £k ,0 < t < T, 



anc 



dW t = /J W s ds, where W t = £ fc e k {t)£ £k ,0 <t<T. 



2.2.1. Wick Product and Skorokhod Integral. 

Definition 3. For £ a , from 3, define the Wick product 



(2-5) ^:=W(^4#)W- 



In particular, taking in (12. 5|) a = kei and /3 = nEi we get 
(2.6) Hki&QHn^) = H k+n &). 
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For a Hilbert space E and arbitrary v = Y2 a v aCa and u = u a ^ a in 
V(b;E), we define their Wick product as a P-generalized real valued ran- 
dom variable given by 



(2.7) vOu = J2J2( u ^ v "-P)z\l py ia G ^'(b; R) 



a /3<a 



Definition 4. Skorokhod integral (Maliavin divergence operator) of v £ 
Li(X>'(b; [0, T], Y)) is a generalized random variable (element of D' (h;H)) 
such that 

fT 



5(v)= [ v(s)dW s = Y j 5(v) a Z 
Jo 



with 

cT 



(2.8) S(v) a = Vv 7 "^ / {v a (k){t),e k (t)) dt. 

k Jo 

and a(k) is given by \2.2\ ). 

live Li(D'(b; [0,T],Y)), then 5 t (v) = f*v(s)dW s = 5 (vl [0>t] ) ,0 < t < 
T, is a process in V (b; [0, T],Y). We have 

St(v) a = ^2Vak {v a (k)(s),ek(s)) Y ds. 
k Jo 

Since W t = Y.k e k{t)£e k , it follows by (J27F]) that 

v a (k) (t),e k (t)) Yy/ak^a , 

a k 

and 

S(v)= [ v t 0W t dt, S t (v) = [ v(s)0W s ds,0 <t<T. 
Jo Jo 

Remark 3. Skorokhod integral is an extension of the ltd integral; \2. 9\) 
below motivates the definition of the Skorokhod integral. 

If u(t) = ^2 a u a (t)£ a is F^-adapted Y- valued such that 

fT 

E / \u(t)\ 2 H dt < oo, 
Jo 

then v = J J u(t)dW(t) = Z*vat 

a is square integrable. By Ito formula for 
the product of JjJ u(s)dW(s) and stochastic exponent pt{z) from Remark [IJ 



Of course, this statement is well known. However, the proof given here is short and 
straightforward . 
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we obtain 

Evp(z) = n[ u(t)dW(t) PT (z)= [ B[p t (z)(u(t),e z (t)) Y ]dt 
Jo Jo 



L 



T 

V[p{z){u{t),e z {t)) Y ]dt,z(L Z. 



So, 

d^Evp(z) ^ f T 9l Q W 



dz 

k 



/*-* Q\^)\ 

Z^ a k ^ Q^)(Vp(z)u(t),e k (t)) Y ]dt, 



(2.9) v a = (V^.)~ l dl '^f (z) \ z=0 = J2[ V^( u a(k)(t),e k (t)) Y ]dt. 

Z k ^° 

Comparing (|2.9p and (|2.8p , we see that Ito and Skorokhod integrals are equal 
in this case. 

3. Wick product Navier-Stokes Equation 
For T > r > 0, let us consider the following Navier-Stokes equation: 
(3.1) d t u(t,x) = di (a ij (t, x) dju (t, x)) + b i (t,x)d i u(t,x) 

- u k (t, x) <>d k u(t, x) + VP (t, x) + f (t, x) 
[<T% x)0iU (t, x)+g (t, x) - VP {t, x)]0W u 
u (r, x) = w (x) , div u = 0. 

The unknowns in the equation (|3.1|) are the functions u = ( u )i<; <( j > P, P- 
It is assumed that a 1 - 7 , 6 l ,f = f/ 1 ) , are measurable deterministic functions 
on [0, oo) x R d , and the matrix (a 1 ^) is symmetric. Let us assume also that 
a ' \ S = (<?*) be y-valued measurable deterministic functions on [0, oo) x R d . 
Let w be a random initial velocity field. 

In addition, we will need the following assumptions. 

Al. For all t > 0, x E R d ,A G R d , 

K\X\ 2 > a i \t,x)X l X j > 5\X\ 2 , 

where K, 5 are fixed strictly positive constants. 
A2. For all t>0,x, 

m&x\d a a ij (t,x)\ + max(|d a bUt, x)\ + \d a aHt,x)\ Y ) < K. 

\a\<2 \a\<l 

A3. The functions f (t, x) and g(i, x) are measurable deterministic, p > d, 
and for all t > 0, 

^[|f(r)|^ + ||g(r)||^r<oo 

l=2,p 

(recall |f(r)| x l7 \\g{r)\\ 1 1 are H£(R ) and HIp(R d , Y)-norms respectively). 
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We will seek a solution to (13.111 in the form 



u(i)=^u a (i)£ Q G P'(b; [0,T],M2),p > 2. 

a 

In this case, denoting by V(v) the solenoidal projection of the vector field 
v, we can rewrite (|3.ip in the following equivalent form: 



(3.2) 8 t u (t) = P[di (a* (t) dju (i)) + b\t)dMt) 

- u k (t) d k u(t) + f (t)] + Pia'WdiU (t) + g (t)]0W t , 
u(r) = w,div u(t) = 0,t G [r,T]. 

If ^ = E a »7a£a with % 6 B*, then (see Remark [2} P(n) = J] a P(r? a )Ca- 

We start our analysis of equation (|3.2jl by introducing the definition of a 
solution in the "weak sense". 

Definition 5. We say that a generalized D-process u(t) = X^a u a(*)£<* £ 
CP'(b; [r,T],H*) is V-W k solution of equation fSDSj} 

in [r, T] , if the equality 

(3.3) u(t)=w+ f V[-u i {s)()d l vi{s) + di{a^{s)djVi{s)) 

J r 

+b\s)d i u(s) + f(s)]ds 

V[a k (s)d k u(s) + g(s)]OW s ds 



holds in P(b; M^- 2 (R d )) /or every r <t<T. If an P-B*- solution in [r,T] 

is also P-H* -solution in [r,T], we call it T>-M k nlj' -solution in [r,T]. 
In the future, we simply say D-solution if there is no risk of confusion. 



Remark 4. L Assume /ioZd, p > 2, w = ^ w a £ Q G P'(b,B*). 

Applying Remark [1| and definition of the Wick product we see that u(i) = 
J^ a u a (£)£ Q G CP'([r, T], b, Hp) is an P-Hp solution in [r,T] if and only if 
for each a, u a G C([0,T], H*) and for t G [r, T] t/ie following equality holds 



2 . 



(3.4) u a (t) = w Q + I V{8 t (a ij (a) ^u Q (a)) + ^(s^u^a) 

- S v (!) ( s ) ^w*) + f ( s ) la =° 

7<« V V// 

+ X] V®ki{ at ( S )> e k(s))y diU a{k) (s) 
k 

+ (g(s),e fc (a)) Y l| a |=i]ds. 
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2. If a = 0, the zero term u a (t, x) = Uo(t, x) of an V-Mp solution in [r, T] 
satisfies Navier-Stokes equation: 

(3.5) u (t) = w + f P[di {a lj (s) ^u (*)) + ^(3)8^) 

J r 

- u k (s) d k u (s) + f (s)]ds. 

For the remaining components we have to solve Stokes equations. For 
\a\ > 1, we can rewrite as 

(3.6) u a (t) = w a + j V[di (o« (s) djU a (*)) + F a (s) 

+ [b\s) - u {s)]diu a {s) - u k a (s) d k u {s)]ds, 

with 



(3.7) F a (s)= Yl jf^VjUyOOW*) 

7<a,|a[-l>|7|>l » ^ 7 ' 

+ V^K^O)' e k(s)) Y diU a{k) (s) + (g (s) , e k {s)) Y l| a | = i]. 

5m ce /or |q| > 1, E£ a = 0, i/ie equation \3. 2\) (or 113.3]) ) can be regarded 
as a random perturbation of the deterministic Navier-Stokes equation h3. 5\) . 

Lemma 1. Let A1-A3 hold, |wo|2,p + | w o|2,2 < 00. Then there is T\ > 
and a unique u G C([0,Ti),H| D 11^) wluwj JX5|) m [0,Ti). 

Proof. According to Theorem 3 in [23], there is Ti > and a unique uo € 
C ([0, Ti) , Hp n H|) such that for each t < Ti 

sup |u (s)|' 1 f+ f |<9 2 u (r)|^r < 00,/ = 2,p, 

0<s<t Jr 

and ()3.5p holds in HT" 1 ,/ = 2, p. By Sobolev embedding theorem, for all 
t<T 1} 

sup |uo(s, x)\ + / sup |Vuo(s, x)\ p ds < 00. 

x,s<t JO % 

and there is a constant C such that for all s G [0, Ti), 

|uo(s)e?fcUo(s)|i, p < C|uo(s)|i,p|9 fe uo(s)| 1)P . 
So, for each i < Ti, 



ng( S )5 fc u ( S )|^<C / K(*)|?^|0*uoMIV* 

j 

< Csup|u (s)|i iP / |9 fc uo(s)|f ds < 00. 

Also, 

ft 

\uo(s)d k u (s)\j 2 ds < 00. 
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Indeed, 

t rt 

\u%(s)d k u (s)\lds < sup \u%(s,x)\ 2 / |Vu (s)||o!s < oo, 

s<t,x JO 



and also 



with 



jT |V (u^( S )d k u (s)) = h+h 

= I |Vu (s)|tds+ / \4(s)d k Vu (s)\ 2 2 ds} 
Jo Jo 

h< sup|Vu (s,x)| 2 |Vu (s)|2ds 

Jo x 

< sup|Vu (s)H / |u (s)|l 

s<t Jo 



p < oo, 



and 



h < sup |u (s,rr)| 2 / \d 2 u (s)\lds 

s<t,x Jo 



< OO. 



By Proposition [5] in Appendix II, u G C ([0,T],H 2 n H 2 .) for every T < 
Ti and ([33]) holds in Lj,Z = 2, p. □ 



Now we fix an arbitrary T < T\ (T\ comes from Lemma [T]) and prove the 
existence and uniqueness of P-solutions to ()3. 2f) in [r,T],r < T. 

Lemma 2. Assume that A1-A3 hold, b G B,w = J2 a w »£a G P'(b;M 2 n 

T/ien /or eac/i r < T < T\ there is a unique P-H 2 n Hg-soZufa'on u(t) = 
Ea«o(*K«eCX>'(b;[0,ri,H?nHi) of in [r,T]. 

Equivalently, for each a, u Q G C ([0, T], H 2 n H 2 ) and /io/d in 

L,,Z = 2,p. 

Proof. According to Remark 01 it suffice to prove the existence and unique- 
ness of a solution for the deterministic system (|3.4p . For a = 0, the existence 
and uniqueness of a solution to (|3.5p follows from Lemma [TJ We proceed by 
induction. Assume there are unique u a G C([0, T], Hp n H 2 ,), |a| < n, such 
that (|3.4p holds in L/,Z = 2, p. By Sobolev embedding theorem, it implies 
that 

(3.8) sup |u a (s, x)\ + / sup |<9u a (s, x)\ p ds < oo, 

x,r<s<T Jr % 

if | a | < n. Then for |a| = n + 1, the equation (|3.6p for u a is Stokes and it 
is readily checked (see ()3.6p ) that 



y |F a (s)|' M ds < oo,Z = 2,p. 
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According to Proposition El there is a unique u tt G C ([0, T], H 2 , n H 2 ,) so 
that ([32D holds in Lj, I = 2, p. ' □ 

Because of the uniqueness, the P-solution has a restarting property. More 
specifically, the following statement holds. 

Corollary 1. Assume that A1-A3 hold, beB,w = E Q ^ ^'(bjE^n 

Let u r ' w (t) 6e i/te uragiie VM 2 p n H 2 , solution to / TO]) m [r,T],T < 2i, 
starting at w. Zei r < r' < t < T. Then 

u r ' w (t) = u r '- u ( r ')(t). 
Proof. Indeed for = u r ' w (t), and r < r' < t < T, we have for £ £ [r, T] 

u(i) = u(r')+ (a ij '(s)^u( S )) +6 i (s)a i u(s) 

-u fc ( S )05 fe u( S ) + f(s) 
+ [^(s^u (s) + g (s) - VP ( S )]0^ s }<is, 
and the statement follows by Lemma [2J □ 
3.1. Rescaling and approximation of the generalized solution. To 

begin with, we will derive more precise estimates for V-M. 2 n H 2 , solutions 
of equation f|3.2|) . One could hardly expect that the P-H 2 n solution of 

unbiased Navier-Stokes equation has finite variance, i.e. ^2 a |u Q (i)| 2 < oo. 
However, in this subsection we will show that the solution could be obtained 
as the limit of square integrable solutions of the equations rescaled in a 
special way using a second quantization operator (see [29j, and Appendix I, 
5-1). 

Fix b = {e n } Gi3 and define an unbounded operator 

Aej- = 2kei~, k > 1. 

Obviously, the projective limit of the domains H„ C H = L<i ([0,T],y) of 
A n with the norm 

I \y\ |h» = \\A n y\ | = (j2 (2^) 2n vl] , y = E v^u e n n , 

\ k / k 

is a nuclear space denoted J\f = Mih). For n G N, let H_ n be the completion 
of H with respect to the norm 

IM|H_„ = \\A^y\\ = (£ {2k)- 2n yl) 1 l\y = Y,yke k G H_„. 

k k 

The inductive limit A/*' = AA'(b) of HL n is the dual of N . For a Banach space 
i£ and p G [0,1], let S P (E) = S p (h;E) be the space of all 77 = J2 a a a£,a £ 
£>'(b,£) such that 

\H\s P , q = H\ P , q = (^(|a|!r(2N) 2 ^|a Q ||) 1 /2 < «, f or every g > 0, 
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where 

oo oo 

(2N) 2qa = Y[(2k) 2qak = 2 2g|Q| Yl k 2qa * 

k=l k=l 

(see [8]). We consider S P (E) = S p (h, E) with a family of seminorms ||f?||5 p;9 . 
Similarly, for p E [0, 1], let S- P {E) be the space of all rj E T>'(h; E) such that 
IMIs_ p _ 9 = ll^ll-p,-? < 00 f° r some q > 0. It is dual of S- P (E) if E is 
Hilbert. If E 1 = R, Kondratiev test function space (AO = 5i(R) and 
Kondratiev distribution space (AT) -1 = <S_i(R) (see [8], p. 39, [6]). 

We will show that the solution found in Lemma [2] belongs to the Kon- 
dratiev space S_i(M2 n H 2 ,) = (M)' 1 (H 2 n H|). 

Proposition 1. Lei A1-A3 hold, sup s fc |efc(s)|y < oo and w E H 2 n H 2 be 

deterministic. Assume that 

u(t) = J2^m a ^ CV'(b;[0,T],M 2 p nm 2 ) 

a 

solve \3.2i) in hi, I = 2, p. Denote 

L a = sup |u a (t)| 2 , P + sup |u a (i)|2,2, a E J. 

t<T t<T 



Then there is a constant B$ such that 
L a < V a 

where K = 1 + supj L £i , and 



!q a i_ 1 W4 Q| - 1 irH,H>2, 



L-iai-i = rn — 7 ii i I' |a| - 2 

11 a — 1 \ a — 1 



are i/ie Catalan numbers (see e.g. [30], [9] j. Moreover, there is a number 
q > 1 so that 

sup u i < ^ i ii < °°> 



i.e., i/ie solution u(t) belongs to the Kondratiev space of generalized random 
functions S_i(IH 2 n H|) = (AA)^ 1 (8) (H 2 n H 2 ,). 

Proo/. For |a| > 1, u a E C ([0,T],H 2 n H|) , are solutions to the Stokes 
equations ()3. 6|) : 



u« (t) = w Q + y V{di (a* (s) d jUa (s)) + F Q (s) 

+ [&*(«) - 40)]diU a (s) - u* (s) d fc u (s)]ds, 
with F a (s) defined by (|3.7p . Finally, by Proposition [5] in Appendix II, 

L a <CY,(f W^Wjdsf' 1 . 



l=2,p 



16 
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Since 

_ 7 (s) S&u T (s)|i ; ; < C|u a _ 7 (s) |2,p|Vu 7 (s)|i 5 i, 

it follows that 



|F a (s)|l,J <C[ J( a )^- 7 ^ 7 



7<a,|a| — 1>|7|>1 
+ Y W0V^A*(fc) + !|a|=l|g(s)|l,i]- 



For \a\ > 2 



L a <C[ ^ yf j Lg-yLj + l^Q ^ ^/OfcL a ( fc )]. 

7<a,l<| 7 |<|a|-l V V// k 

So, there is a constant Bq so that for \a\ = n > 2,L a = (a!)~ 1//2 L a . 
L=.. we have 



\7<a,l<|7|<|«| — 1 k 

Denoting L a = L a if \a\ > 1, L a = 1 + L a if |a| = 1, we get 

^« < -Bo L Q _ 7 L 7 

-y <Q;, 1 < | 'y | < | ct | — 1 

and by [9] for |a| > 2 

i 



with 



So, 



K = l + C[Lo + ^(rilBWIlW 7 ']- 

l=2,P J ° 



11 < a!q 2 ahl ( l " l )(2Nr J B 2(|Q| - 1) ^ 2 H 



and 



^L 2 

a 



* < C | 2 Qhl ^ l ^(2Nrr J B 2(|Q| - 1) K 2 l Q l. 
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Therefore with r = (n),ri = (2i)~ 2q ,q > 1, 



|a|=n |«|=n 



For large n, the Catalan numbers 

AP—l 

Cn-1 



V^F(n - I) 3 / 2 
and there is a number q > 1 such that 

n=0 |o|=n 



Therefore, the solution u(t) belongs to Kondratiev's space *S'_i(lHIpnIHl|). □ 

Remark 5. // Let u(t) is the solution \3. 2\) with a deterministic w 6 M. 2 n 
Hi, and q > 1 is the number in Proposition (TJ i/ien f/ie action of u(t) = 
^ a u a(t)^ a can be extended from V(h) to 



'1,9 



W = = ^OaCa € 5i(R) : \\f]\\s ltq < Ooj 



<7.s 



(u(t),r/) = ^u a (t)a a ,r? = a Q £ a G <Si, g (R). 

a a 

Note that the stochastic exponent p{z) = p(e z ) (see Remark\T\) belongs to 
5i ;(? (R) in (a) provided 

\A q e z \ 2 Y = 5> fc 2*fc*) 2 < 1. 

k 

For e > define a self-adjoint positive operator D £ on such that 
D e eu = 2~ ek ek and a sequence of positive numbers k £j11 = e~ ee " . Set 
C e = J2^=o K e,nDf n . It is a second quantization operator in the Fock space 
U = Y.n^Un = H® n (see Appendix I, 5.1). Clearly, 

(3.9) C £ e a = K Et{a] Df n e a = K £j]a] (2~ eN ) a e a , 

where 

oo 

(2- £N )° = J\2~ £ka K 

k=l 



18 R. MIKULEVICIUS AND B. L. ROZOVSKII 

Proposition 2. Assume A1-A3 hold and sup s k \ek{s)\y < oo. Let u(i) = 
E a u a0)6* G CT>'([0,T],b,H 2 nH|) 6e a generalized V-UlnB%-solution of 
equation 113. 1\) in [0, T] wii/i a deterministic w G Hp n H| and 

oo 

U £ (i) = C e u(t) = ^ Ua(t)(2- £N )^ Q , 

n=0 |a|=re 

where C £ u is rescaling based on the second quantization operator C £ (see 
Appendix I, 5.1). 

Then u e (i) is Mj^-valued square integrable process satisfying the equation 
(3.10) d t u £ (t) = V{8i [a ij (t) dju £ (t)) + ^(t)^u £ (t) 

- c £ (c £ - l u k £ (mta^dku^t)) + f (*) 

+ CeitfWC^diUe (t) + g (tyOiC^Wf)]}, 

u £ (0) = w, div u £ = 0. 

Moreover, u £ (t) G 5i(H 2 D Hi) = (AT) 1 ® (H 2 n H|), t G [0,T],u(t) = 
C £ u £ (t) and 

sup||u £ (i) -u(t)|| 5 _ ^0 

t<T 

as e — > 0, where q is a number in Proposition [TJ 

Proof. LetL Q = sup t < r |u Q (t)| 2iP +sup t < r |u a (t)| 2 , 2 - Since u £>Q (t) = K e ^ (2~ £N )° u Q (t), 

L £jQ , = SUp \u £>a (i)\2,p + SUp |u £iQ (i)| 2 ,2 = K e ,| a | (2~ £N )° P Q . 
t<T t<T 

Since for each q' > 0, there is a constant C(e, qr', g) independent of a so that 
(|a|!) 2 (2N) 2 ^4 H (2~^) a 
<^(2fc) 2 ^ 2- 2efe )H e -2-H ( | a | !) 2 

< C(e,q',q) < oo, 
it follows by Proposition [1] that 

IIUeWHl, . < C(S, q', q) £ ^L!!!^ < oo. 
So, u £ (t) G Si(H 2 nH|),t G [0,T]. In particular, 

E|u £ (t)|| >2 = K,a(t)\h < ^ M- 

a 

Therefore u e (t) is H|-valued square integrable process and (|3.10p follows by 
Remark [TU] in Appendix I, 5.1. Obviously, u(i) = (C £ ) u £ (f). Since 

u £ (t) - u(i) = J> - (2- £N ) a ]u Q (t,x)^, 
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it follows that 



2 (2N) 



-2qaf2 



u e (t)-u(i)||L 1 ,_ ! <El 1 -^H ( 2 " N ) 



a 



-> 



Q: 



as e — > by Lebesgue's dominated convergence theorem uniformly in t. □ 

Remark 6. The solution u m Proposition [JJ depends on a fixed uniformly 
bounded basis b in H. belongs to the Kondratiev space 



constructed using a Gelfand triple M C H = L2 ([0, T],y) C A/ 7 which 
depends on b. 

In [7], a class Q^ 1 C (M)" 1 of regular generalized functions was introduced 
that does not depend on a fixed Gelfand triple or a basis in H. Unfortunately, 
the estimates in Proposition [JJ (because of the factor N _<?a = ITfcA; _ ' 3afe / ) do 
rao£ imply that u is a regular generalized function of class Q -1 . Also, the 
space 



with some q > 1, to which the solution in Proposition [JJ belongs, cannot be 
embedded into any space with weights depending only on \a\ (for example, 
into the spaces, like Q~ l ~ q in [7], that do not dependent on a fixed basis or 
Gelfand triple in H): inf| Q |_ n (N)~ qa = 0. 

4. Markov property and independence of basis 

In this Section we will show that a generalized T>-Mp n Hg-solution of 
equation (|3.2|) has the following properties: it is adapted with respect to the 
filtration (J 7 ^) generated by the Wiener process Wf, it is independent of 
the choice of the basis b, and it is a generalized Markov process. 

4.1. Equivalent characterization of D-generalized processes. A more 
convenient characterization of P-solution to (|1.2p (see Definition [5]) is based 
on another (equivalent) description of T>(h). It allows to introduce the notion 
of an adapted solution and extend it from V(h) to a space of test functions 
that is independent of b G B. 

4.1.1. Equivalent description of test function space. Often it is convenient 
to use the exponents p(z),z G Z, defined in Remark CD to describe the test 
function space T>(h). 

To each multi-index a of length n we relate a set K a whose elements are 
positive integers ki,i = l,...,n, such that each k is represented there by 
ctfc-copies. An ordered n-tuple K a = {k\, . . . , k n } with k\ < k2 < • • • < k n 
characterizes the locations and the values of the non-zero components of a. 
For example, k\ is the index of the first non-zero element of a, followed by 
max (0, ttfej — 1) of entries with the same value (see [21]). 



5_i(m2 n ml) = (My 1 ® (m 2 n e|) 



S-i,- g (Bj n Ml) = {r/G5_i(HjnHi) : \\v\\s^. q < 00 



} 
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For an orthonormal basis {ek,k > 1} in L2 ([0,T],Y) and a £ I with 
K a = {ki, . . . , k n }, we denote 

E a =Yl e K w ® • • • ® e fc.(„) » « G J > 
where G™ is a permutation group of {1, . . . , n}. The set 

(4.1) \e a = agjl 

[ Va!|a|! J 

is a CONS for the symmetric part H®" of H® n . 
For \a\ = n, 

(4.2) i a = ^\d\\W{e a ), 
where 

(4.3) W(e a )= / .../ e tt (ai ) ...,a n )dW n ... ( iW Sn . 

If a = e k , then W(e £k ) = W (e fc ) = / Q T e fe (t)dW t . 
According to f)2.4p . 



/II! 

p(z)=p(e 2 )=p(z,b) = ^ ^ ^yei W (e a ) 

n=0 \a\=n 

= ^ ^ 7^' 

n=0 | a |=n V 

Denote 

Pn(z) =Pn(e z ) = V ~ =£a 

1 va! 

|a|=n 

(4.4) =/ / •••/ e«(« 1 )...e a (s„)dW 81 ...dW r anJ n>2, 

Jo Jo Jo 

Pi(z) = Pi(e z ) = / e z (si)dW sl ,p (z) = p (e z ) = 1. 
j 

Lemma 3. For b Ze£ V =V(b) be the linear space of random variables 
that consists of all finite linear combinations ofp n (z),z £ Z, n > 0. T/ien 
V(b) = £>(b) fin particular, £ Q G V(b)j. 

Proo/. Obviously, V = V(b) C P = V(h). For a <E 1, denote «(«) = 
max{/c : a& 7^ 0}. Fix iV, n and a = («&) 6 I such that |a| = N,k(o) = n. 
Consider a finite dimensional Hilbert space 

G = \ ^2 Va ^ a : Va G R \ 

I \a\=N,K(a)<n I 
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with inner product 

(Y^ U <*£a, 2 V '^ G = Y VaV '*- 
a a a 

Let G be a vector subspace of G generated by pw(z), z = (z%, . . . , z n , 0, . . .) G 
Z. It is enough to show that G = G. Indeed, the subspace G is finite- 
dimensional and obviously closed. Assume there is a vector v a ^ a G G 
which is orthogonal to G. So, for all z = (zi, . . . , z n , 0, . . .) G Z, 



z 

{y2v a ^ a ,p N (z)) G = Vw a7 = = o 



which implies that all v a = 0. Therefore G = G.This completes the proof. 

□ 

Due to Lemma O we can characterize convergence in T>' = T>'(h) by test 
functions of the form p m (z). Indeed, for z G Z, v G P, and m > 0, we have 

(4.5) (Pm(z),v)= V W a ^=. 

i | va! 

Therefore we have the following necessary and sufficient condition: 



Corollary 2. A sequence v n — > u in 2?' i/ and on/y if for all z £ Z and all 
m > 

(p m (z) , u n ) -> (p m (z) , u ) . 

4.1.2. Action of a Skorokhod integral onpM(z)- Consider v(t) = Y2 a v a(t)^ a £ 
T>' (b; [0, T], K) such that for all a,/c, J* T | (u Q (s), efc(s))y |ds < oo. Recall 
that the Skorokhod integral assigns to such v a generalized random process 



[ v(s)dW s = 6 (vl m ) = V 6 t (v) a Z a ,0 <t<T, 
Jo 

St(v) a = Y,V®k {v a (k)(s),e k (s)) Y ds. 
k ^° 

Remark 7. For pm{ z ) G T>(h),z £ Z, M > 1, it is easy to show that 
(pm{z) ,S t (v)} = [ ( (pm-i 0) ,v(s)} ,e z (s))y ds. 



4.1.3. Action of a Wick product on pm(z)- Recall that for a Hilbert space 
E and arbitrary v = v a ^ a and u = J2 a u a ^ a in V(b, E), we define 
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In particular, 

W/3)! 



(3\a\ 



The following statement holds. 

Lemma 4. For a Hilbert space E, arbitrary elements v = v a ^ a and 
u = Ua^a from V(h, E), and z S Z, M > 0, 

(p M (z),v0u) = (pk(z),v) (p L {z),u) ; 

K+L=M 

In particular, 

{l,vOu) = ((l,v),{l,u)) E 
(the generalized expected value of v()u is the product of expected values). 

Proof. According to (14. 5p . 

(p M (z),v0u) = EE}^^=^^ 



\a\=M P<a 



E E -m u ^ E ~m u ^ 

K K+L=M \/3\=K VH ' | 7 |=L V ' 

((p{z),v) , (p{z),u)) 



E 



□ 



4.1.4. An equivalent characterization of the solution. Now, we will charac- 
terize the solution of equation (13. lh by its action on test functions pm(z), z £ 
Z, M > 0. The following statement holds. 

Remark 8. Assume A 1- A3 hold, w = E a w «C« € P'(b, H^nH|), divw = 
and u(f) = £ a u a (i)£ a G CP'(b; [r,T],ll2 n M|). Then u(t) is V-U 2 p n H| 
solution of \3. 1\) in [r, T] if and only if for all z £ Z and M > 0, 

u^(i) = <u(t) >PA r(z)> = ^ u " ,/,i " 



\a\=M 

is continuous in M. 2 n Hi and i/ie following equality holds in hi, I = 2,p, 
(4.6) u A/ >* (i) = w M - z + fv[di (a ij («) d jU M ' z (s)) 

J r 

+ 6 i (s)5 i u M ' 2 (s)- J2 u k > K » z (s)d k u L > z (s) 

K+L=M 

+ l M=0 f ( s ) + {a l (s), e z (s)) Y d t u M -^ z (s) 
+ lM=i(gO) ,e z (s)) Y ]ds, 

where M > 0, w M,z (x) = (w(x),pm(z)) = Yl\a\=M w a (x)z Q /v / aI, and 
u _1 ' z (t,x) = 0. 1/ M > 1, equation ^4.6\ ) is Stokes equation; if M = 0, 
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it is Navier-Stokes equation. Indeed, we obtain fl^.ffi ) by multiplying both 
sides of j3.4\ ) by z a /\/~a\ and adding. 

4.2. Adapted and independent of basis generalized processes. Let 

Loo ([0, T] ,Y) be the space of measurable Y- valued bounded functions on 
[0, T). For h G Lqo ([0, T],Y),M > 0, we denote 

PM,t(h)= / ... / h( Sl )...h(s M )dW Sl ...dW SM ,0<t<T. 
Jo Jo Jo 

By ([HI), p M (z) = pu(e z ) = PM,r(e z ),z £ Z. 

Lemma 5. (i) If {m k ,k > 1} is a CONS in L2 (0,T) , and {£k,k > 1} is a 
CONS in Y ,h € Lqo ([0, T] ,Y), then for each n,n' > 1, there is z £ Z such 
that 

n' n 

K,n'(t) = ^2^2 / (h(s),£ k )mi(s)ds£ k mi(t) = e z (t), 
i=l k=i ^° 

< t < T. Obviously, PM,T(h n ,ri) G ^(b), b = {e^ = mi k £j k , k>l}, 

h n ,n> ^ h in L 2 ([0,T],y), 
PM,T(h n>n i) ->■ pu,T(h) in L 2 (fi,P), 

as n, n' — >■ 00. 

(ii) Assume (rrii) is trigonometric basis or unconditional L p ([0, T]) -basis 
(for example, Haar basis, see |19j ). /i G Loo([0, T], Y). Then there is a 
sequence z(n) £ Z such that e 2 ( n ) —> h in L p ([0,T],Y) for all p > 2, as 
n — > 00. 

Proof. We prove the second part of the statement. Let 

n 

h n{s) = '^2(h(s),£ k ) Y £k,n > 1. 

k=l 

Then |/i n (s)|y < |/i(s)|y and for all p > 2, 

f \h n (s) - h{s)\ p Y ds -> 
Jo 

as n — > 00. If (mi) is trigonometric basis or unconditional L p ([0, T]) -basis 
(for example, Haar basis), then for each n 




K,n'{s) - h n (s)\^ds -> 



as n' -> 00. So, there is a subsequence Z n such that 

/ \h n ,i n (s) - h(s)\ p Y ds ^ 
Jo 

as n — > 00, and (ii) follows according to part (i) of this remark. □ 
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Let T be the space of all linear combinations of PM,T(h), h G L oo ([0, T],Y), 
M > 0. Obviously, UbG^f (b) C 7", and T does not depend on any particular 
be B. 

We say that h n -> /i in ^([O.TJ.F) if /i n -> /i in L p ([0,T],Y) for all 
P > 2. 

Definition 6. Given a Banach space E, let T' = T'(E) be the space of all 
linear E-valued functions v on T such that v(pM,T{h n )) — > v(pM,T(h)) for 
all M > if h n — >• h in L oo ([0, T], Y). We say v G T'{E) is a generalized 
E-valued r.v. 

Denote T(b) = T'(b,E) the set of all v G V(b,E) such that for each 
h G L oo ([0,T],y) and any sequence e z r n \ — >• h in Loo([0, T], Y), the limit 
linin^oo (pm,t i z ( n )) ) w) exists in E for all M > 0, and does not depend on 
a particular sequence z{n) such that e z ( n ) — > h in L oo ([0, T], Y"). We define 

(PM,T (h) , v) = lim (p M (z(n)) , u) . 

n— >oo 

Lemma 6. Zei beB. XTten 

fa; r'(b) cr'(s). 

(b) For v G T'(E), there are a a £ E such that the restriction 
v\v(b) = J2 a ^M £T'(b,E). 

a 

Proof, (a) Let v G 7~'(b), h n — > h in L oo ([0, T], Y). For each n there is a 
sequence 2fc(n) such that e. Zk ( n ) ~^ h n in Loo[0, T] and f(zfc(n)) — > v(h n ) as 
— > oo. Therefore there is a subsequence 2fc n (n) such that 

l e ^„(n) ~ ^lin([0,T],Y) + l^(e^„(n)) - v{h n )\ < 1/n. 
Then for each p > 2, n > p, 

\tz kn (n) ~ h\ Lp ( [0 ,T],Y) 

< \ e z kn (n) ~ K\l p ([0,T],Y) + \ h n ~ % p ([0,T],Y) 

<Tp » \e Zkn{n) - /i„U n ([ ,T],y) + \K~ Hl p {[o,T\,Y) 

So ' e ^„W h m L °° (t°' T ]' y ) and 
\v(h n ) - v(h)\ 

< \v{h n ) - v(e Zkn(n) )\ + \v(e Zkn(n) ) - v{h)\ -> 

as n — > oo. 

(b) Let W Gr , ( J B),b = {e n },e a = ^(b),7V> 1, 

Pn{z) =p N (z,h) =p N {e z ) = V — ^=f£q- 

M=7V 

By Lemma O £ Q G 7" and by linearity 



>(pn(z)) = V -P=v(£ a 
, i «, v«! 



|a|=JV 
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Denoting v a = v(£ a ) set v = Y2 a v a £,a G £>'(b). Obviously, 
v{pn{z)) = ^2 -1=^^) = v{p N (z)) 

\a\=N 

and (b) holds. □ 

Let T(b; [r,T]) = T'(b; [r,T],E) be the space of all v G P'(b; [r,T],E) 
such that v(t) G T'(b,E),r <t<T. 

Definition 7. Given a Banach space E, let T' = T'([r,T], E) be the space 
of all T' (E)-valued functions v(t) on [r,T]. We say v G T'([r,T],E) is a 
generalized E-valued stochastic process. We denote CT'([r,T], E) the set of 
all continuous u G T'([r, T], E). 

The following obvious consequence of Lemma [6] holds. 

Corollary 3. Let b G B, r < T. Then 

(a) T'{b;[r,T],E) Cf([r, T], E). 

(b) For v G T'([r, T], E), there are E-valued functions a a (t),0 < t < T, 
such that the restriction 

v\v(h) = J> Q (i)£ Q (b) G T'(b;[r,T],E). 

a 

Now, we introduce the notion of an adapted generalized process. 

Definition 8. (a) We say v G T'(E) is J-^ -measurable if 

(pM,T(h),v) = (p M ,t (h),v) 

for all /i G ^ ([0, T] , Y) , M > 0. 

(b) We say v G T'([0, T],E) is ¥ w -adapted if v(t) is -measurable for 
each t. 

Example 2. Let Wt be a cylindrical Wiener process in a Hilbert space Y 
and Wt = -^Wt- Then (see Example^ as well) Wt and Wt are generalized 
Y -valued adapted stochastic processes. For any h G L oo ([0,T] , Y), 

(Wt, PM ,T(h)) = f h{s)ds = (Wt, PM ,t(h)) , 



[Wt,PM,T(h)) = h (t) = (Wt,PM,t( h )) 
ifM = I, and (W u pf{h)) = (w t ,pf(h)) = 0, < s < T, otherwise. 



4.3. Independence of basis and Markov property of the solution. 

Remark [8] suggests the following definition of a generalized solution to (|3.2p . 



Definition 9. Given w G T'(H^ n M|),divw = 0,T > r, a generalized 
process u G CT'([r, T],B^f] H|) is called H% n M^-solution of equation [3J]) 
in [r,T], if for each h G Loo([0, T],Y),M > 0, the function u M ' h (t,x) = 
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(pM,T(h), u(t, x)) is an Ml n M^-valued continuous functions satisfying in 
hi, I = 2,p, Stokes (M > 1) or Navier-Stokes (M = 0) equation 

(4.7) u M ' h (t) = w M ' h + J' V[di (a ij (a) o^u^ (a)) 

+ 6 l ( S )9 l u M '' l ( S )- ^ (a) Qfcii^OO 

A'+-L=M 

+ l M=0 f (s) + (^(s), fc( S )) y ^-^(s) 
+ l M =i(g (s) ,h(s)) Y ]ds, 

where M > 0, w M ' h (x) = (p M ,T(h),w(x)) , and u~ 1 ' /l (t,ac) = 0. 

Obviously, a generalized solution is a D-solution. Now we are in a position 
to prove the main result. 

Theorem 2. Assume that A1-A3 hold, w G T'(H£ D Hf). T/ien /or eac/i 
r < T < Ti tfiere is a unique Hp 1 n U^-solution u G CT'([r, T],B% n Hf) 
o/ ( fff.l?]) in [0, T]. Moreover, if w is -measurable, then the solution is 
Y w = [Jt ) t> -adapted and it extends all T>(h)- solutions, b G B. 

Before proceeding with the proof of Theorem[2]we shall prove the following 
auxiliary statement. 

Lemma 7. Let A 1- A3 hold, w G V (Mj n H|) , div w = 0,/t G ([0,T], Y"). 
TTien i/ie infinite system 7| j wi/t Mv M,h = (w,PM,T(h)) has a unique solu- 
tion v M ' h G C ([r,r],H^ nH|) ,M > 0. Moreover, /or eac/i iV > 1 tfiere is 
a constant C independent of h so that 

(4.8) i?^<C{ R kRl+ Y^ N,K \^ + Rn -^ f\Ks)\ l ds) 1 ' 1 

K+L=N, l=2,p ^ r 

K,L<N-1 

+ l N =i( [ T \h{ S )fds)W{ f MsW^dsf' 21 }}, 

J r J r 

where R N = sup r < s < T [\v N > h {t)\ 2 , p + |v^(i)| 2 , 2 ] . 

Proof. If M = 0, then (|4.7|) is Navier-Stokes equation (j3 . 5[) and, by Lemma 
[T] there is a unique Hp n H^-valued continuous solution v° = uo- We proceed 
by induction. Assume there are unique v A/ = v M ' h G C([0,T],H^ n H|) 
solving the system for M = 0, . . . , N — 1. Consider the equation for v^: 

(4.9) v N (t) = w N ' h + I P[di (a« (s) fyv* (s)) + ^(s^v^s) 

Jr 

- (s) d k v N (s) - ^- fc (s) d fc v°(s) + ffc(s)](fo 
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with 

fft(s) = - Yl yK ' k ^ d ^ L ( s ) + lAf=i(g (a) , fc(a))y 

K+L=iV, 
K,L<N-1 

+ {a i (s),h(s)) Y d t v N - 1 (s) 

Since 

/ |f^(s)|^ds<oo,Z = 2,p, 

the existence and uniqueness follows by Proposition Also, by Proposition 
[5j for each N > 1 there is a constant C independent of h such that, denoting 
R% = su Ps < T [\v M < h (s)\ 2;P + W M > h (s)\ 2 ,2] ■ Since for I = 2,p, 

\^(s)h,i <c\ \v K 00 Mv L (*)k P + |M*)|y I^Wki 

iC,£<iV-l 

+ |fc(*)M|g(*)lli.i, 

the inequality (|4.8|) follows by Cauchy-Scwarz inequality. □ 

4.3.1. Proof of Theorem^ Fix T < T\ and choose a special CONS b = 
{e„} G £> such that for each h G Loo([0, T], Y) there is a sequence z(n) G Z 
(see Lemma[5|) for which e x t n \ — > h in L p ([0, T],Y) , for all p > 2, as n — )• oo 
(for example, taking in L2([0,T]) a trigonometric basis or unconditional 
L p ([0, T])-basis (Haar basis), see [19J). By Lemma[6]with £ a = £ a (b), 

w i©(b) = Wa ^°- 

a 

According to Lemma there is a unique V-Ml n M|-solution 

u(t, x) = Y u « (t, x)& G CD' (b; [r, T] , Hp n H|) 

a 

of (|3.2p in [0, T]. The coefficients u a (t,x) satisfy (|3.4p and, by Remark [U 
fOjJ) holds for all M > 0, z G Z. Fix h G Loo([0, T],Y) and consider an 
arbitrary e z t n \ — > h in L p ([0, T],Y) , for all p > 2, as n — )• oo. 
Then 

w M,z(n) = (w, PAf (z(n))) 

= w Q (t,2;)z(n) a /\/d GM^nM^, 

|o|=M 

and 

u^W(^) = (u(M), PM (z(»))) 

= u Q (t,x)z(n) a /\/^! GC([0,T],H^ne|), 

M > 0, is the unique solution to the system (|4.7|) corresponding to /i = e z / n ) 
and w M ' /l = w M ' z ( n ). Recall, u _1 ' z(n ) = and u 0,z(n ) coincides with the 



2S 
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solution of Navier-Stokes equation uo in Lemma [TJ By Lemma El there 
is a unique v M ' h £ C ([r, T], n H|) , M > 0, solving ([471 with w M ' h = 
(w,p M (h)) G Hp n B%. We have v~ 1,?i (t) = 0, and v°' h coincides with the 
solution of Navier-Stokes equation Uo in Lemma HJ By Lemma [7J (see (|4.8|) ). 
for every M > 1 there is a constant Co independent of n such that 



(4.10) 



sup [l^'^Wk^+lv^Wkd < C . 



r<t<T 



l=2,p 



For V 



M 



r M,h 



u M,z(n)^M > 1, the following equation holds (M > 1) : 



(t) = w M - w M > 2 W + / (o lJ (a) a,-V« (s)) 

+ [6*( s ) - u (a)]0 i V* f (a) + V^(a)5iUo(a) + G n ( S )]ds, 
where G n (s) = G*(s) + G£(s) with 

GiW = - E [^(^u^s) + (a) 9 fc V^( s )] 



K+L=M, 
K,L>1 



and 



G*(M) = lM=i(gOO, M*))y + (^( S ),/ i ( S )) y 5 i Vf- 1 ' /l ( s ) 



By Proposition[S]in Appendix II, for L\ = sup s < T | V^(f)| 2iP +| V n {t)\i,i 
with M > 1, 



where A r . 



L n M <C[A n +^(f IGnWijds) 1 / 1 }, 

l=2,p Jr 

J2i=2, P l wM - ™ M ' z(n) \2,i- We estimate 
|Gi( a )| lf , < CC £ Lfc, 



1<K<M-1 



|G=(ft,s)|i,,<C7[lM=i|/in(a)|y|g(a)|i I i 

+ ^5tf_i + C7oM*)l], 
where h n = h — e z r n \. So, for each M > 1 there is a constant independent 
of n such that 

L « f < C {A n + Yl L K + E [ / T 
l<if<M-l «=2,p ^° 

+ (^ T |^( S )| 2 ^ S )^ ( || g ( S )|f^ S )^] } . 

Starting with M = 0, Lq = for all n, it follows by induction that 
L n M = sup E|v^*( a )-u M '*W( S )| 2 ,,->0 



r<s<T 



i=2,p 
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as n — > oo, M > 1. 

Since h G Loo([0, T],y) is arbitrary, 



u(t) = ^u a (t)e a g cr'Cbj^ri^nei) c CT'([r, r], Hp n Ml) 



is the unique solution to (J3JJ). Obviously, 

v M ^(t,z) = (p M ,T(h),u(t,x)} ,h g Loo([o,r],y) 

satisfies (|4.7p . Since for each CONS b' = (e^,) G any linear combination 
of e' k belongs to L oo ([0, T], y), the generalized solution u(i) extends any 
P-solution. 

Now we will prove that the unique generalized Hp n H^-solution of (13.11) in 
[0,T] is F w -adapted. We fix t* G (0,T),r < t* and consider a special basis 
b G B with rhiit) in L2((0,T)) so that each fhi is supported either in [0,t*] 
or in [t*,? 1 ] and such that for each h G Loo([0, T],Y) there is a sequence 
z(iV) G Z (see Lemma [5]) for which — > h in L p ([0,T],y), for all 
p > 2, as iV — > oo (for example, (m^) is a combination of two trigonometric 
or unconditional L p ([0, T])-basis (Haar basis) on (0,i*) and (t*,T)). Let 
Cq = la(b), a G J, the corresponding orthonormal basis in L2(F^). Let 

u(t)=^u a (i)f a GCr / (b;[0,T] ) H2nH|) 

a 

be the unique solution to (|3.ip constructed using the representation 
w| c(S) = w a £ Q G T'(b; H^ n Mi). 

a 

So, u a G C([0,T],H^nH|) satisfy {33} in IL/, Z = 2,p, with w Q = w Q ,a G 
J. Let J' ={a G J : a has a non zero component corresponding to fhk 
whose support is in (t*,T)}. Since w is J-^-measurable, w a = if a G J'. 
Indeed, if a G J', there are q G R, Zj S = l...,n, so that 

n 

£a = ^QPt(^)- 
i=l 

In \ n 
Wo = ( W,^CjPr(^) ) = ^Cj (w,p T (zi)) 

\ i=l I i=l 

n I n \ 

= y^Cj (W,p r (^)) = / W,y^CjP r (^) \ = 0, 
i=l \ i=l / 

because, by (|4.3|) . 

n n 

i=i i=i 



Then 
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We claim that, similarly, for t E (0, t*), u a (t) = if a E J'. Indeed if 
a E J' and \a\ = 1, we have w a= and in the equation ()3.6f) for u a we have 
the input function F a (t) = for t E (0, t*). Therefore the unique solution 
= if t E (0, t*). Then we simply apply induction on |a| = n and use 
(|3.6p (note that if \a\ = n + 1, then a = a + for some k and without any 
loss of generality we can assume that a E J'). As a result, 

u(t) = E u Q (t)e Q = E * e [r,t*]. 

Obviously, £ a are J-^-measurable for a ^ J'. Also, for any z E Z, M > 
l,t<t*, 

(p M {z),u(t)} = E Ua{t)—r= = {pM,t*{z),u(t)} 



= (pM(e z l(Q,t*)),u(i)) 

(note that e z = X^fc z fc e fe) e zl(o,**) = Ylk^G z k e k, where G is the set of all k 
such that fhj k in = rhj k lj k has its support in (t*,T)). The statement of 
Theorem [2] is proved. 

The solution above has the restarting property as well. By the same 
arguments as in Corollary [1] we have 

Corollary 4. Let w E T'(H|nEl|) and A1-A3 hold. Let u r - w (t) be H£n]H|- 
solution to h3. l\) in [r,T],T < T%, and r < r' < t < T. Then 

(4.11) u r ' w (t) = u r '- u ( r ')(t). 

Corollary 5. (Markov Property) Assume that the assumptions of Corollary 
\Qhold true and, in addition, w is —measurable, then vT'™ (t) is (J~t)t>r- 
adapted. This together with {1^.11 ) can be interpreted as Markov property. 

Let choose a uniformly bounded basis b = {e n } ( sup fc s |efc(s)|y < oo) 
and rescale the solution in Theorem [2] 

u ll>(b) = y^Ug(t)^q, 

a 

using the second quantization C £ operator in f|3.9[) in the Fock space H = 
Yun U n {Kn = H®", see Appendix I, 5.1). Recall, 



n=0 



where D £ ek = 2 ek ek and k £iU = e ee ". According to Proposition [21 

oo 

(4.12) rf(t) = C £ u(t) = E E ^,n(2- £N ) Q u a (t)e Q 

n=0 \a\=n 

is Hg-valued continuous. The following statement holds. 
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Proposition 3. Let b = {e n } be uniformly bounded and A1-A3 hold. Let 
u G T'([0,T],e^ nHl) 6e £/ie solution to i TO) in [0,T] u«f/i tfie determin- 
istic initial w G Hp n H|. T/ien i/ie rescaled M^-valued square integrable 
continuous process u £ (t) (defined by ^.Pty , see Proposition^ is adapted 
and Markov (in a standard, rather than generalized, sense). 

Proof. For each M > 1, z G -Z, 

Ep TjM (z)u E (i) = (p T ,M(^),U £ (t)) = K £iM (PT,Af(^),u(t)) 
= Ke,M (pt,M(^ £ ),u(t)) = (p ttM {z),U £ (t)} 

= B PtiM (z)u £ (t). 

Since u e (i) is square integrable, E [u £ (i)|.7 : ^n = u e (i). So u £ (t) is adapted 
in a standard sense. 

For any < s < t, by Corollary U u(t) = u s ' u ( s )(f). Therefore, 

u e (t) = C £ u(t) = C £ (u s ^ s \t)) = (u E ) s ^^(t) 

and the standard Markov property follows. □ 

Acknowledgement 1. We are very grateful to S. Kaligotla and S. Lototsky 
for useful discussion. 

5. Appendix I. Wiener chaos 
In the first part of Appendix we present some facts of white noise analysis. 

5.1. Rescaling of Wiener chaos by second quantization operator. 

Consider a generalized random variable u = Y^ a u aCa = Ua\/\c*\W(e a ) G 
V'(h), b = {e k , k>l}eB, where 

W(e a ) = W* n (e a )= / .../ e a ( Sl ,...,s n )dW Sl ...dW Sn , 
Jo Jo Jo 

if \a\ = n and {e Q , a G J} defined by (|4.ip is a CONS of the symmetric part 
U n = H® 71 of H® n (recall H = L 2 (0,T]) x Y). We can interpret 

oo 

a n=0\a\=n 

as a result of the noise W acting on an element of the Fock space: W(u) = u 
with 

oo 

u = y^u a ^/\a\\e a = u a Vn\e a 

a ra=0|o|=n 



oo oo 



n=0 n=0 n=0 
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Here H®° = R and the norm in the Fock space H is defined as 

oo oo ||-(n)||2 

|2 _ || V^ f .(n)|,2 _ST^ W U \Wn 



it 

n! 

n=0 n=0 



Obviously, E[W({i) 2 ] = Let A = (A n ) n>0 be a self-adjoint positive 

operator in H such that A n e a = X (a) e a , where |a| = n and A (a), a £ J, 
are positive numbers. 

Remark 9. T/te operator A in the Fock space H can be used to rescale a 
generalized r.v. For 



u 

a 

we define Au = u A G ^'(b) by 



Au = u A = W(Au) = ^u a ^/\a\\W{Ae a ) 

a 

= ^u a \{a)^/\a\\W \e a ) = ^2u a X(a)£ a . 

a a 

Definition 10. Since A(a) > 0, we can define 

A- l u = u A ^ = X n « A («) _1 ^- 

a 

Example 3. 1. (second quantization operator in space-time) Consider a 
self-adjoint positive operator B in H such that Be k = Xk^k- The second 
quantization operator A =T(B) = (5® n ) in % is defined as 

Ae a = T(B)e a = B® n e a = X a e a , \a\ = n, 

where A = (A&) and X a = IlfeA^*. We have 

T(B)u = ^2u a X a ^ a . 

a 

2. (second quantization operator in space) Consider a self-adjoint positive 
operator B on Y such that the sequence of its eigenvectors (£ p ) p>1 (B£ p = 
Xplp, X p > 0) is a CONS in Y. Let b = {e^, k > 1} , where ek(s) = mi k (s)£j k . 
We extend B to H by 

Be k = B(m ik £ jk ) = m ik B£ jk = X jk e k 

and rescale in space-time using A = (B® n ). For u = Yla u a^a we have 

T(B)u = Y J UaX a ^ a , 

a 

where X a = II fe A? fc . 

Jk 
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3. Consider a self-adjoint positive operator B on H such that Bet = X k e k 
and a sequence of positive numbers q n . Let A = Ylri=o QnB® n . Then Ae a = 
q n B® n e a = q n \ a e a , \a\ = n.In this case, 

oo 

Au=^2q n ^ u aX°£ a . 
n=0 \a\=n 

For the Wick product we have the following obvious statement. 

Remark 10. Assume A = (A n ) is a self-adjoint positive operator on % such 
that Ae a = A n e a = X(a)e a , \a\ = n and X(a),a G /, are positive numbers, 
u, v G X>'(b). Then, denoting Au = u A ,Av = v A , we have 

A(uOv) = A (A-WOA-W) =Y, C »^, 

a 

where 



EX(a) a A a ' 

<Q A(/3)A(a-/3) M ^^\/ /?!(«-/?)!• 

In particular, if T(B) is the second quantization operator in space-time, 
then 

Xja) _ X a _ 
A(/3)A(a-/3) ~ XPX a -P ~ 
andT(B)(uOv) = u B 0v B . 

2. For the Skorokhod stochastic integral, we have 

A (5(u)) = A J Q u s 0W s ds = J A (A^ufOA^Wty ds 
A (A^ufOWs) ds; 



for the coefficients 

A(a) 



(A6(u)) a =J2[ (^ {k) (t),X(e k )e k (t)) Y dty/c^ 



X(a(k))X(e k ) 



" T X(a) 



X(a(k)) 

(A(u(t)0W t )) a = J2(v^ k) (t),e k (t))y^^-^ a . 

k 

5.2. Product, Wick product, and Malliavin derivatives. Consider 
Hilbert space H = L 2 ([0,T\,Y), CONS b = {e k , k > 1} G B, cylindrical 
Brownian motion Wt, and Cameron-Martin basis {£ a } Qg j introduced in 
Section [2.1i The Malliavin derivative T> (see e.g. [28]) is defined in T>(b) 
as follows (it assigns to £ a an element of T>(h; H)): 

(5.1) P(^) = J2V^^- £ k e k = Yl ^2 V^k e k^a- 

k>l a ft=a+e k 



a • 



X(a) 
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By induction, 



(5.2) V^) =J2[ E WW ^5^«p I 



\p\=n 



where 



u p = Ek u ...,k n E £fel +...+ £fcn = p efci ® ■ ■ • ® e fc „ € H® n . 



Proposition 4. Lei and £ re 6e elements of the Cameron- Martin basis. 
Then, with probability 1, 

(5.3) = ^ . 

n=0 



Proof. It is a standard fact (see e.g. [18J) that 



P/ VP/ V K — P 



Let us rewrite this expression as follows: 



^ V^!(/3 + 7 )! , 

p+7=K,p+/3=6 

where the summation goes over all triples (p, /3, 7) G Jx J x J such that 
p + 7 = K,p + P = 6. Changing variables (1-to-l mapping that assigns to 
(p, (3, 7) the vector (p,/3,a) with a > /?) of summation by p = p, f3 = (3, 
7 + (3 = a, we get 

(5.4) = E ^Sv^ 



a 



pW(a-/3)! 

p+«-^=K,p+/3=6» 

E E E E w-nW=« /3)i e 



Q 



By definition of the Wick product and (15. 2p . and taking into account that 
is the number of different orthogonal unit vectors in u p G H® n , we 
arrive at 

(5.5) V n Ze<>V n ti K 



-1^2^ nl 2^ ip+ a -0=«ip+P=o Kma _ /3)l Z«- 

a fi<a \p\= n 

Comparing (|5.4p with (|5.5p . we get f|5.3j) . □ 
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Remark 11. Proposition \4\ implies that £$£ K = Ce'O'^.k + Y^^<e+K C 7^'r- ^ n 
other words, ^eO^fc = (,9+k is the highest stochastic order component of the 
Wiener chaos expansion of £q£ k . 

By linearity, the statement of the Proposition could be extended to 



(5.6) XY = J2 



(V n X)0(V n Y) 



n=0 

where X and Y are finite linear combinations of elements of Cameron- 
Martin basis. If X and Y finite second moments, relation 15. 6\) could be 
derived from the former case by passing to the limit in L\.By linearity, the 
statement of Proposition could be extended to 

xy _ (V n X)0(V n Y) 
n\ 

n=0 

where X and Y are finite linear combinations of elements of Cameron- 
Martin basis. 

5.3. Derivation of unbiased Navier-Stokes equation. To simplify dis- 
cussion, we will consider a velocity field which depends only on one stan- 
dard Gaussian random variable rj ~ iV (0, 1) , rather than a trajectory of 
the Wiener process Wt- An interested reader would have little difficulties 
extending the arguments below to the setting with Wiener process. 
Consider a velocity field 



u (*! x ) = ^2 u ™ x ) • 



n=0 

Note that in our setting the Cameron-Martin expansion (see Theorem [I]) is 
indexed by integers rather than multi-indexes. Assume that for every n, u n 
is analytic in x in that it could be written as 

u n (t,x) = c„ i7 (t) x^ 1 . 

7GN d 

Let Z = {Z\, . . . , Zd) be a ^-measurable. Then by substituting Z into u 
we get 



(5.7) 



n \ 7 / 



Now, let us introduce the Wick-powers of Z : Z^ := zf 11 ()...() Z 1 / , 7 = 
( 7 i,...,7d) GN d . 

Next we will replace the standard algebra in (|5.7p by the Wick algebra: 
u£ (t,Z) := J> n , 7 (t)Z^ 
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Consider now the following random field 

U (t,Z) := J>£ (t,Z) 0Cn(v) 

n>0 

Remark 12. Note that Wick algebra on nonrandom elements reduces to the 
standard deterministic algebra. 

Let X t = (Xj , . . . , Xf ) be a solution of the following dynamic equation 

X t = u° (t,X t ). 

Then by the Wick chain rule 

X t = ^-u (t, X t ) = dtvfi (i, X t ) + Vu° (t, x t ) ox t , 
dt 

= dtvP (t, X t ) + u° (t, X t ) V0u<> (t, x t ) . 

If F = F(t,x) is an acting force, this yields (Wick) Euler equation 

dtu (t, x) = -u° (t, x) V0u° {t, x)+F (t, x) 

If there is no randomness, due to Remark 1121 this equation reduces to the 
standard Euler equation: 

d t u (t, x) = — u (t, x) Vu (t, x) + F (t, x) . 

Now, by taking F = Au — VP, where P stands for pressure, we get the 
unbiased Navier-Stokes equation 

<9<u° (i, x) = Au - u° (t, x) VOu (t, x) — VP + F (t, x) . 

6. Appendix II. Stokes equation 

Consider a deterministic Stokes equation for u =(n')i<;<^, and scalar 
functions P, 

d t u(t,x) = di (a' lj (t,x) dju(t,x)) +b' l (t,x)di\i(t,x) 
+G(t,x)u(t,x) +f (t,x) + VP(t,x),divu(t) = 0, 
u (r, x) = w(x),x 6 H d , r < t <T. 

equivalently, 

(6.1) d t u (t) = S[di (o« (i) fyu (i)) + 6 4 (t)a,u(t) + G(t)u(t) + f (t)] 
u(r) = w,t G [r, T], 

where 5 is the solenoidal projection of the vector fields, 

a(t,x) = {aVfaxfl^^Mt,*) = {b l (t,x))^ t<d ,G(t,x) = {g^{t,x)) x ^ 

are measurable bounded functions. The matrix a is symmetric and positive. 
We will need the following assumption. 
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B. For all t > 0, x, A eR d , K\X\ 2 > a ij (t,x)X i X j > 5|A| 2 ,where K,S are 
fixed strictly positive constants. Also, for all (t, x) G [r, T] x R , 

max \d%a(t,x)\ + max \d£b(t, x)\ + sup \G(t, x)\ < K. 

t,x,\a\<2 t ,x,|a|<l t,x 

Definition 11. A function u G C([r, T], W s p (R d )) is an M s p - solution of l6l\) 
if the equality 

(6.2) u(t) = w + S[d i (a i ^s)d j u(s))+b i (s)d i u(s) + f(s)]ds, 
ZioWs in EIp _2 (R d ) /or ever?/ i G [r, T]. 

Proposition 5. Let p > d, assumption B /ioW ; w G M p n H|, 

/ T |f Wli.jda + / T |VG(a)|jcfe < oo, Z = 2,p. 

TTien i/iere is a unique MpHM^-valued continuos solution to \6.1\) . Moreover, 
there is a constant C independent of f , w, u so that 

(6.3) sup [|u(i)| 2 , 2 + |u(i)| 2jP ] 

r<t<T 

< c (|w| 2>2 + | W | 2iP + \m\i P dsf P + mihds) 1 ^ . 

Proof. Let w G 1^(1 Hf. By Proposition 4.7 and Corollary 4.6 in [22] 
(applied for s = 0), there is a unique Hp n H 2 -valued continuous solution u 
of HSU, and 



/ l u (s)| 2 ,«^ < oo,/ = 2,p. 



Moreover, 



(6.4) sup|u( S )|^<C[|w| 2i ,+ [ T \i(s)\ l l ds],l = 2,p, 

s<T Jr 

Consider Stokes equation 

(6.5) £(*) = w a + J t V[d i (a^(s)d^{s))+F(s)]ds 
div £(t)=0,t€[r,T\, 

where 

P(s) = b\s)diu(s) + G(s)u(s) + f (a). 
It is readily checked (using Sobolev embedding theorem) that 

(6-6) |F(a)|i,, < C[K(\u(s)\ 2il + | VG(s)|i|u(s)|i iP 

+\i(s)\ 1 , l ],l = 2,p. 

and 

/ \F(s)\l l ds<oo,l = 2,p. 

Jr 
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By Corollary 4.6 and Proposition 4.7 in [22J (applied to (|6.5f) with s = 1), 
there is a unique Hi n Hp-valued continuous solution of (|6.5j) £ = u (by 
uniqueness) such that 

[T 

/ |u(s)ls iZ ds < oo,/ = 2,p. 

Let a be a multiindex such that |a| < 2. Then u a = d a u is LpHL^-valued 
continuous and satisfies the equation 

d t u a (t) = S{d a [d i (a i '(t)d J u(t))+F(t)], 

u Q (0) = d a w. 

Differentiating the product, we obtain 

(6.7) d a di{a i] {t)d 3 vi{t)) = d l {a lj {t)d J vL a {t))+di^o i {t) 
with 

(6.8) \T) a (t)\ l <C\u{t)\ 2>h l = 2,p. 

By Lemma 3 in [25], yi, a {t) = |u a (t)|[,Z = 2,p, is differentiable: 

yi, a {t) = yiA r ) + / h iA s ) ds ' 

J r 

with 

h l:a (s) = l{(\u a (s)t 2 u a (s),d a F(s)]) 1:l 

a^(s)d i (\u a (s)\ l -\ k a (s))d j u k a (s)dx 



Notice d a F(s) G H-i^ and, by our assumptions, there is a constant C so 
that for all s E [r, T] 

|9 a F( S )|_ liZ <C|F( S )| M ,/ = 2,p. 
We have h^ a { s ) = lhj a (s) + lhf a (s), where 



= - / a i ^s)d i (\u a (s)\ l - 2 u k a (s))d J u k a (s)dx. 

Then 

^»<-<5 1 |u Q ( S )|'- 2 |Vu( S )| 2 dx, 
and for each e > there is a constant C e such that 
|fc?«(*)l<e / \u a (s)t 2 \X7u a (s)\ 2 dx 

+ C7 e |[|u a ( S )r 2 (|VF( S )| 2 + |D Q ( S )| 2 ) 
+ |u Q (s)|'- 1 |F( S )|]dx J 
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So, we obtain that 

2/i(t) = ^ = ^( r ) + / hi(s)ds 

\a\<2 Jr 

with 

/»,(*) = E fyaW^faOO + /«(*))> 

|a|<2 

where /;(s) = |F(s)|^ + |D a (s)|j. Therefore, by (|6Tg]) . JETS]) . 

< + |fWli,j + |VG( S )|||u( S )| z 1)p . 

By Gronwall's inequality, 

T 



sup y p (s) < C[y p (r) + / |f( S )|? )P ^] 



and 



(6.9) ^up r |u(t)|^ <c(K jP + ^ |f( S )|? >p d S ) , 

where C is independent of w and f . Similarly, by Gronwall's inequality 
sup y 2 (s) <C[y 2 (r) + sup |u(s)|? iP f \VG(s)\ 2 2 ds + f |f (s)|? >2 ds], 

r<s<T s<T Jr ir 



and ( see ([6 
(6.10) 

sup |u(i)|i )2 <cf|w| 2 + |w|| p + (/ \f( s )\P ds )l + [ \m\hds 

r<t<T \ Jr Jr 

where C is independent of w, f and u. Combining (|6.9p and (|6.10p we have 
(|6.3p with C is independent of w, f and u. 

Given w G Hp n Hi, there is a sequence w n G Hp 1 n H 2 so that w n — > w in 
Hp n H 2 . For every n there is a unique Hp n H 2 -valued continuous solution 
u n of (|6.2p with the initial condition u(r) = w n . By (|6.3p 



V] sup \u n (t) - u m (t)\ 2 ,i < C } |w„ - w m | 2ji -)■ 

l=2,p — — '=2,p 

as n, m — )• 00. There is a continuous Hp n H 2 -valued u(t) such that 

V] sup |u n (i) - u(i)| 2 ,j -> 

l=2,p — — 

as n — >• 00. Obviously, u is Hp n H 2 -valued continuous solution of (|6.2p with 
initial condition u(r) = w and (j6.3[) holds. □ 
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